Power spectrum of many impurities in a d-wave superconductor 
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Recently the structure of the measured local density of states power spectrum of a small area 
of the Bi2Sr2CaCu2 0g (BSCCO) surface has been interpreted in terms of peaks at an "octet" 
of scattering wave vectors determined assuming weak, noninterfering scattering centers. Using 
analytical arguments and numerical solutions of the Bogoliubov-de Gennes equations, we discuss 
how the interference between many impurities in a d-wave superconductor alters this scenario. We 
propose that the peaks observed in the power spectrum are not the features identified in the simpler 
analyses, but rather "background" structures which disperse along with the octet vectors. We 
further consider how our results constrain the form of the actual disorder potential found in this 
material. 



In the past few years, high-resolution scanning tunnel- 
ing microscopy (STM) experiments on the cuprate su- 
perconductor BSCCO QSHmilllillla, IISII3 

have obtained local information on electronic structure 
for the first time. The first great success of this tech- 
nique was the observation of resonant defect states at 
low temperatures in the superconducting state P, |3, 13 , 
confirming early proposals that such states should be 
reflected in the local density of states (LDOS) of d- 
wave superconductors [T^ IT^. Subsequent experiments 
revealed the existence of nanoscale inhomogeneities|3, 
0, 0, which are currently the subject of debate, be- 
ing attributed either to interaction-driven effects such as 
stripe-formation0,|3 or to Friedel oscillations of weakly- 
interacting quasiparticles0, 1^. At the heart of this de- 
bate lies the important question of whether, and if so 
in what ranges of doping, conventional BCS-like theories 
describe the ground state and low energy excitations of 
BSCCO. Recently it has been pointed out that, in inho- 
mogeneous systems, the Fourier transform of the LDOS 
(FTDOS) contains information not only about the dis- 
order potential, but about the kinematics of the associ- 
ated pure system. That this must be so at some level 
is clear from the original one-impurity problem solved 
by Friedel a charge inserted in an electron gas gives 
rise to LDOS oscillations which vary at large distances 
as ~ cos 2kFr/r^, so the wavelength of the LDOS "rip- 
ples" caused by a single impurity gives the Fermi wave 
vector directly. A somewhat more sophisticated version 
of this arg ument,^10, 15, ]J,[Il[l3 still assuming scat- 
tering from a single or few impurities and noninteracting 
quasiparticles, suggested that the peaks in the FTDOS 
are due to scattering of quasiparticles by a weak disorder 
potential. In this case favored momentum transfers cor- 
respond to vectors q connecting two tips of quasiparticle 
constant energy contours which maximize the joint den- 
sity of states (JDOS). This interpretation has been ap- 
plied to recent experiments which claim to map out a 
Fermi surface in agreement with angle-resolved photoe- 
mission (ARPES). In contrast to this, Howald et al.|l|i| 
identify nondispersing features in the FTDOS from their 



experiments, which they suggest are indicative of static- 
stripe formation along the Cu-0 bond directions. 

In this Letter, we report on numerical studies of models 
of disordered superconducting BSCCO. Previous analy- 
ses of a single weak impuritv fisl ll^ agreed qualitatively 
with the simple JDOS analysis of Ref. and succes- 
fuUy predicted the dispersion of the peak positions in 
the FTDOS. It has been since shown jl9j that the sharp 
1-impurity peaks survive in the many-impurity case if 
the potential is sufficiently weak. However, the ability to 
correctly predict the peak positions (which depend only 
on the quasiparticle band structure, and not the disorder 
potential) does not imply that the low-energy excitations 
are understood. A full microscopic understanding of the 
superconducting state requires an accurate description of 
the disorder potential in the BSCCO Cu02 planes. We 
will show that there are both qualitative and quantita- 
tive problems describing experiments which stem from 
the inadequacy of the weak-impurity model, and that 
the FTDOS spectrum of a reahstic model comprising a 
dilute concentration of unitary scatterers, together with 
a smooth disorder potential, is required to fit experiment 
and has features not present in the simple JDOS analy- 
sis. Understanding these discrepancies will be crucial to 
extracting reliable information about the clean system, 
as well as the nature of the disorder potential. 

Many Impurities. In this section, we discuss the dif- 
ferences between the single and many-impurity FTDOS 
for weak pointlike scatterers, and compare these spectra 
with experiment. We first imagine a random distribu- 
tion of Nj pointlike impurities located at sites Ri with 
i = l...Nj. The LDOS can be formally decomposed 
p{r,uj) ~ po{uj)+6p{r,uj) where po is the DOS of the ho- 
mogeneous superconductor, and Sp is the local shift due 
to disorder given exactly by 



6p(r,uj) 
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^Im ^ [G"(r - R.)4g°(R, - r)] ^(^1) 



where the w-dependence is suppressed for clarity, T(uj) is 
the 2Nj X 2Nj many-impurity T- matrix (the factor of two 
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FIG. 1: FTDOS at lo =14 meV for weak potential scat- 
ters (Vb = 0.67ti): (a) for one weak impurity, with a few 
important scattering wavevectors indicated; (b) for 0.15% 
weak scatterers. Cuts through the data of (a) (thick line) and 
(b)(thin line) along the (110) direction and scaled by 1/V Ni 
are plotted vs. in (c), while (d) shows the weak scattering 
response function Im A3(q,aj). Peaks at q = are removed 
for clarity. In all figures, the x and y axes are aligned with 
the Cu-O bonds. 



arises from spin), G°(r,a;) is the bare electron Green's 
function (^^(r, w) = G°(k, w) cxp(ik • r), " refers to 
matrices in the Nambu-spinor formaUsm and [. . are 
spinor indices. The T-matrix is expressed in terms of the 
1-impurity T-matrix i, = [1 - ViG^{r = 0)]-^V, by 



Ni 



]G'°(Ri — Rm)Tmj, (2) 



m=l 



where the impurity potential at Ri is Vi = VqTs, and fi 
are the PauU matrices. For pure potential scatterers, the 
T-matrix can be decomposed into two Nambu compo- 
nents t = toe^'^°To + t^e^'^^Ts with to and ^3 real, and (po 
and 03 the scattering phase shifts. 

Throughout this work, we adogt the ARPES-derivcd 
tight-binding model of Normanj20| for the band struc- 
ture, and assume nearest-neighbor d-wave pairing with 
order-parameter Ak ~ Ao(cosfc^ — cos ky) and Aq — 
0.16^1 = 24 mcV where ti = 150 meV is the nearest- 
neighbor hopping matrix element. We take, as represen- 
tative, the weak scattering potential Vq = 0.67U ■ We 
ignore the effects of tunneling matrix elements pa. l23l |. 
and assume that the STM probe measures the LDOS 
p(r, Lu) directly (note that r need not correspond to a 
site R of the crystal lattice). The Fourier transform is 



then p(q, w) = X^reixL^ *'''''p(r, w), where L x L is a 
square set of positions at which measurements are 
made, and q = 27r(m,n)/L are vectors in the associated 
reciprocal lattice. 

The FTDOS is shown at w = 14 meV in Fig. [H for (a) 
a single impurityQ (b) a collection of Ni = 21 (0.15%) 
weak scatterers on a 120 x 120 lattice. Using the notation 
of Ref . [Tll | , a few important peak positions which are pre- 
dicted by the JDOS analysis are also shown in Fig.^a). 
To understand the relationship between Figs, ^a) and 
(b), we consider Eq. to leading order in f^: 

6p{q,uj) =i V t„(q)Im[e'^°A„(q,c^)] (3) 



a=0,3 



where <„(q) = ta e"''*-^ 
Ek [G°(k,w)f„GO(k 



and Aq, (q, uj) = 

q, a;) . Eq. ||2Jl affords a clear 

separation between degrees of freedom associated with 
the disorder potential and those of the pure system, 
which is not possible if scattering processes of higher or- 
der in ti are important. Since ta{ci) consists of a sum 
of Nj random phases, it becomes a random function of 
q in the first Brillouin zone as Nj — » 00. By contrast, 
Aa{q,Lu) is the response function of the clean system, is 
independent of the disorder potential, and determines the 
peak locations and widths in the FTDOS. As is evident 
in Fig.^c), JDOS peaks are not broadened or shifted by 
disorder; note in particular that the thin line represent- 
ing the many-impurity case shows a sharp q7 peak which 
is visible only as one or two bright pixels in Fig. ^b). 

In the limit of weak potentials, Eq. (PJ reduces to the 
result of Capriotti et al[T^. 



Sp{q,uj) ~ -V{q)Im A3(q,w)/7r, 



(4) 



which is also valid for finite range V{r). The ef- 
fect of an extended potential can be understood 
with an example: if the single-impurity poten- 
tial is Vi3(r) = Wo exp(— r^/2r§), then V^(q) = 
27rrQUo cxp(— q^rQ/2) E- exp(iq • Ri). Thus, for sufh- 
ciently weak scatterers, the q-space structure of (5jo(q, uj) 
is determined primarily by the band structure, but has 
an envelope which suppresses the FTDOS near the Bril- 
louin zone edges, and is noisy because of randomness in 
the impurity positions. We remark that since A3(q, a;) 
can be calculated from knowledge of the electronic struc- 
ture of the pure system, Eq. Q is, in principle, a useful 
tool for determining the scattering potential directly from 
experimental measurements of the LDOS. However, we 
caution the reader that Eq. applies only to unphysi- 
cally weak potentials. Figure^d) shows that even for the 
weak potential Vq = 0.67ii, not all peaks in the FTDOS 
(notably the central peak around q = 0) are reproduced 
by Eq. Q. 

While Eqs. Q and appear to suggest that the dis- 
persing 1-impurity peaks are also present for many impu- 
rities, we can show from Eq. Q that in the disorder aver- 
age (|Mq,^)l') Ni and {\6 piq, u;)\^) - {\S piq, iu)]^ ^ 
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Ni{Ni — 1), implying that the noise is actually as large 
as the signal. As a consequence, the weak qi peaks 
present in Fig.^a) are lost in the noise in Fig.^b), and 
the broader background features are relatively enhanced. 
Thus, it appears that the peaks predicted by the JDOS 
analysis may not be robust when many impurities are 
present. This is particularly true when the effects of en- 
ergy resolution Auj arc considered. In Ref. Aw w 2 
meV|27j. and we use a complex energy uj + in our 
calculations which for 7 = O.OlSfi = 2.25 meV yields 
a comparable resolution. We find that even for one im- 
purity, FTDOS peaks are extremely sensitive to 7. All 
peaks are suppressed, but not broadened, as 7 increases 
(this conclusion differs from Ref. 

It is therefore difficult to understand, based on the 
weak impurity analysis, why experiments measure broad 
and well-defined qi and qy peaks of roughly equal weight. 
Furthermore, 1-impurity calculations such as in Fig.^c) 
find in addition to JDOS peaks a highly structured dis- 
persive "background" structure, whereas in current in- 
terpretations of experiments all features are ascribed to 
JDOS peaks. Our calculations suggest that, because of 
noise, it may be easy to confuse dispersing background 
features and JDOS peaks in experiments. Very recently, 
it was pointed out in Ref. that the scattering ampli- 
tude for processes involving Ao(q) is not sharply peaked 
at the "octet" qa, but our analysis shows further that 
the predominant structures in the presence of noise are 
generally considerably shifted from the qa positions, and 
may indeed be associated with off-shell processes. 

Realistic Disorder Models. The weak-scattering point- 
like disorder model is convenient for its simplicity, but 
does not predict FTDOS peak widths consistent with 
experiment. In addition, the relative weights of the fea- 
tures ascribed to the q^, (e.g., the nearly equal weight 
of q7 and qi) are not correctly predicted, particularly 
at low energies. In this section, we explore whether the 
measured FTDOS can be explained by a simple "realis- 
tic" disorder model. Two types of disorder, unitary de- 
fects and intrinsic nanoscale inhomogeneties are known 
to be present in nominally clean BSCCO. In the ex- 
periments of Ref. ^3 1 ^ concentration of w 0.2% uni- 
tary, pointlike native defects were observed, while in 
Refs. 0, IE @j S B 13 nanoscale inhomogeneities were 
observed with a typical size of « 2 nm. In this work, 
we model the native defects as pointlike scatterers with 
an on-site potential Vq = 30ti which produces a lo- 
cal resonance centered at meV (the concentration and 
strength of the unitary defects is thus fixed by experi- 
ment within the potential scattering model). The source 
of the nanoscale inhomogeneity is unknown, and we take 
the simplest ansatz, that it arises from a smooth ran- 
dom potential, probably originating from charge inhomo- 
geneities in the BiO layers. For dcfiniteness, we model the 
smooth potential as V{r) = X^i ^(*) 6xp(— ri/A)/fi and 
fi = [(r - R;)^ + d^]^^^, where Rj -|- zd^ are the defect 
locations, V{i) are the defect potentials and A is a screen- 
ing length. We take a bimodal distribution V{i) — ztV so 
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FIG. 2: Comparison of theory and experiment. Calculated 
LDOS and FTDOS are compared with experimental FTDOS 
(EXPT) at five energies. LDOS are shown for a 40 x 40 sub- 
section of the 120 X 120 lattice. FTDOS are shown in the first 
Brillouin zone, — 7r/a < qx,(ly < Tr/a. Theoretical calcula- 
tions are for a mix of unitary pointlike and smooth potential 
scatterers. Scattering g-vectors are shown at 22 meV and 14 
meV (negative energy vectors are the same). A broad back- 
ground feature in the FTDOS is identified with an arrow at 
-14 meV. 



that the smooth potential represents spatial fluctuations 
of the local potential about a mean which is determined 
by the doping level. Equation (0) applies to pointlike 
impurities only, so we use an implementation of the re- 
cursion method [2^ to solve directly for the local Green's 
function G(r, r,uj + 17), with 7 = 2.25 mcV, of the inho- 
mogcneous system. Our best-fit model consists of 0.2% 
unitary defects, and 8% smooth scattering centers with 
V = 2ti, d= 2a and X = a. 

We compare this mixed-impurity model with the ex- 
perimental data of McElroy et al in Fig.|21 In both exper- 
iment and calculations, the local resonances are plainly 
evident in the LDOS for ^15 mcV and, in our calcu- 
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lations, make the dominant contribution to the FTDOS 
As with weak scatterers, the narrow JDOS peaks are 
swamped by noise, and suppressed by the finite energy 
resolution. However, unhke the weak scattering case, the 
qi and structures at 14 meV are relatively robust 
and remnants of the peaks can be seen. Excess weight 
relative to experiment in the Brillouin zone corners is 
probably due to finite potential range and sub-unit cell 
information obtained in experiment, both of which lead 
to a modulation of Sp{q) which decays at large qp^. 

For larger energies, the LDOS impurity resonances dis- 
appear in experiment but, because the electronic wave- 
functions are excluded from the defect site, they remain 
clearly visible in our calculations. The discrepancy may 
arise because the STM tip height is adjusted to main- 
tain a constant current at a fixed bias voltage (typically 
-150 mcV in Ref. ^3), so that what is actually plot- 
ted is the LDOS relative to the average LDOS at that 
site. We will discuss the effect of tip-height adjustments 
elsewhere p^. and here we simply note that despite dis- 
crepancies, the Fourier transforms of the measured and 
predicted patterns are qualitatively similar. At energies 
\uj\ <; 15 meV, the smooth random potential makes a 
significant contribution to the FTDOS. 

The agreement between the calculated and experimen- 
tal FTDOS is qualitatively good at low energies, although 
our calculations show more asymmetry between positive 
and negative energies, with the negative energy calcu- 
lations fitting better, than in experiments. This is the 
result of the large, probably artificial, asymmetry of the 
model band, in which there is a van Hove singularity at 
w ~ — 50 meV that is not observed in tunneling exper- 
iments. As before, remnants of JDOS peaks (for exam- 
ple, the qi peak at -14 meV or the qs peak at 14 meV) 
are evident in the calculated FTDOS, but most of the 
structure comes from a set of broader dispersing features, 
which has been previously ignored. At -14 meV and -22 
mcV, we see that there are clear octet structures about 
the central peak, but that the dispersing "q7" features 
[the arms along the (±7r, ±7r) directions extending from 
the central peak] extend well beyond qy. In general, the 
background structures disperse qualitatively (growing or 
shrinking with \uj\) as one expects from the JDOS anal- 
ysis. 

Since the q7 peak comes from intranodal scattering, 
it is a direct measure of the k-dependence of the su- 
perconducting gap and scales with ^ I/^a, where va 
is the gap velocity at the nodes. The measured gap is 
not consistent with nearest-neighbor d-wave pairing — 
for example the experimental value at lu = —14 mcV 
q^^^ w O.GTr/afll] is nearly twice the value in our nearest- 
neighbor d-wave model — and the authors of Ref. have 
been forced to introduce a significant subleading cos 69 



harmonic in their fit of the angular dependence of the 
gap on the Fermi surface. The additional harmonic is as- 
sociated with pairing beyond near-neighbors, in contrast 
to ARPES results which suggest the gap is pure near- 
neighbor d-wave at optimal doping. Recognizing that the 
observed feature at roughly twice the true qj is in fact 
the background feature (similar to that seen as a "hump" 
in Fig. 1(c)) found in our calculations may enable one to 
bring the two experiments in closer agreement. 

At w = ±22 meV, the calculated (110) structures, 
which correspond to forward (intranodal) scattering, are 
stronger than the (100) structures. Preliminary numer- 
ical calculations suggest that scattering from order pa- 
rameter fluctuations may contribute to qi-type peaks, 
leading to the interesting question of how the interplay 
of disorder and order parameter fluctuations manifests 
itself in the FTDOS. It is a general feature of our calcu- 
lations at higher energies that the qi peaks are weaker 
than observed in experiments, and a tendency to stripe 
formation along the Cu-0 bonds, as proposed in IS 1 
could possibly enhance the weight near the qi peaks. 
However, stripe formation is difficult to reconcile with 
the dispersal of the qi peakfTollTH. 

Conclusions. Existing analyses of recent STM mea- 
surements of the FTDOS have established the fundamen- 
tal point that the structures in the FTDOS are dispersive 
and therefore likely represent interference of disordered 
quasiparticle wavefunctions. We have argued, using a 
combination of analytical and numerical approaches, that 
previous weak-scattering analyses are inadequate to ex- 
plain the details of the FTDOS, however. In previous 
works it was expected that sharp "octet" peaks would 
broaden into the observed structures due to impurity 
scattering; here we have shown that these peaks are lost 
in the noise created by the interference of many impuri- 
ties, rather than broadened. Our work suggests that the 
dispersing features observed may not correspond directly 
to the predicted sharp peaks at all, but rather to dis- 
persive background structures which have hitherto been 
ignored. We show that the simplest "realistic" disorder 
model of unitary impurities and a smooth random poten- 
tial explains many features of the FTDOS at lower ener- 
gies, but at higher energies the comparison is worse, prob- 
ably refiecting our lack of knowledge of the true smooth 
component of the disorder, including the local order pa- 
rameter fiuctuations neglected here. Further refining of 
these comparisons will be extremely important in un- 
derstanding the origin of the nanoscale inhomogeneities 
observed in STM experiments. 
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